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Abstract 

The propagation of damage in a confined magnetic Ising film, with short range competing magnetic fields 
(h) acting at opposite walls, is studied by means of Monte Carlo simulations. Due to the presence of the fields, 
the film undergoes a wetting transition at a well defined critical temperature Tn,{h). In fact, the competing 
fields causes the occurrence of an interface between magnetic domains of different orientation. For T < Tw{h) 
(T > Tw(h)) such interface is bounded (unbounded) to the walls, while right at Tw(h) the interface is essentially 
located at the center of the film. 

It is found that the spatio-temporal spreading of the damage becomes considerably enhanced by the presence 
of the interface, which act as a "catalyst" of the damage causing an enhancement of the total damaged area. 
The critical points for damage spreading are evaluated by extrapolation to the thermodynamic limit using a 
finite-size scaling approach. Furthermore, the wetting transition effectively shifts the location of the damage 
spreading critical points, as compared with the well known critical temperature of the order-disorder transition 
characteristic of the Ising model. Such a critical points are found to be placed within the non-wet phase. 

PACS numbers: 75.70.-i, 75.30.Kz, 75.10.Hk, 05.10.Ln 
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I. Introduction 



The study of the thermodynamic properties of confined systems has attracted considerable attention over the last 



decades [2 [21 E [1 [5l [6l H [HI El [H [III [H [El [H [111 [IS [13 [H The critical behavior of 



confined fiuids is rather different from the bulk criticality due to the subtle interplay between finite-size and surface 
effects. The interaction of a saturated gas in contact with a wall or a substrate may result in the occurrence of very 
interesting wetting and capillary condensation phenomena, where a macroscopically thick liquid layer condenses at 
the wall, while the bulk fluid may remain in the gas phase [23 [23 [211 [3Ql [HI]. The wetting of solid surfaces by 
a fluid is a phenomena of primary importance in many flelds of practical technological applications (lubrication, 
efficiency of detergents, oil recovery in porous material, stability of paint coatings, interaction of macromolecules 
with interfaces, etc. [271I321I331I34] ). 

Wetting transitions are also observed when a magnetic material is confined between parallel walls where com- 
peting surface magnetic fields act. For example, when an Ising [35l[36] film is confined between two competing walls 
a distance L apart from each other, so that the surface magnetic fields {H) are of the same magnitude but opposite 
direction, it is found that the competing fields cause the emergence of an interface that undergoes a localization- 
delocalization transition. This transition shows up at an L— dependent temperature Tyj{L, H) that is the precursor 
of the true wetting transition temperature Tw{H) of the infinite system [HI HHl [Sl [HZj- It should also be remarked 
that, although the discussion is presented here in terms of a magnetic language, the relevant physical concepts can 
rather straightforwardly be extended to other systems such as fiuids, polymers, and binary mixtures. 

On the other hand, the study and understanding of the propagation of perturbations in magnetic material is also 
a subject of increasing interest. For this purpose, the damage spreading method is a powerful and useful technique 
that has been applied, for example, to Ising systems [HHl [HSl [iQl [43 [121 [IHl [lH [Ml [46l ^\ as well as to spin glasses 
[48l[49j, for a review see e.g. [50j. 

In order to apply the damage spreading method one has to start from an equilibrium configuration of the magnetic 
material at temperature T. Such a configuration is generically called the reference or unperturbed configuration 
S^{T). Subsequently, a perturbed configuration S^{T) is obtained just introducing a small perturbation into 
S^{T). This procedure can be achieved, e.g. by fiipping a small number of spins of the unperturbed configuration. 
The time evolution of the perturbation is of primary interest because, either, it can become healed after some time 
(eventually it may remain finite) or it can propagated over the whole system. The latter scenario may be undesired 
when a thermally stable magnetic material needs to be achieved. A simple measure of the perturbation is the 
"Hamming distance" or damage between the unperturbed and the perturbed configurations |38[ 139). Then the time 
evolution of the perturbation can be followed evaluating the total damage D(t) defined as the fraction of spins with 
different orientations, that is 



where the summation runs over the total number of spins N and the index I, {1 <l < N) is the label that identifies 
the spins of the configurations. 

In our previous study of damage spreading using Ising magnets confined in two-dimensional geometries [46] and 
in absence of external magnetic fields, it has been shown that the presence of interfaces between magnetic domains, 
in the direction perpendicular to the propagation of the damage, causes the spatio-temporal enhancement of the 
spreading. Due to this interesting feature and in order to contribute to the understanding of the role played by the 
interfaces, the aim of this work is to investigate the propagation of the damage at the interface generated by the 
localization-delocalization transition of the Ising magnet confined in presence of competing magnetic fields. In this 
case, and in contrast to the previous work [46], the propagation of the damage occurs in the direction parallel to 
the interface. In order to carry out the study, Monte Carlo simulations of the Ising magnet in confined, thin film, 
geometries are performed. Imposing open boundary conditions and surface magnetic fields along the surfaces of the 
films, the propagation of the damage in both directions, parallel and perpendicular, to the domain interfaces can be 
studied. Furthermore, extrapolations of the results to the thermodynamic limit allow us to determine the critical 
points of damage spreading and locate them in the wetting phase diagram, as evaluated exactly by Abraham [37] . 

The manuscript is organized as follows: in Section II the numerical procedure for the simulation of damage 
spreading in confined geometries is described. Section III is devoted to a brief discussion of the equilibrium con- 
figurations characteristic of confined geometries with competing applied magnetic fields. The results are presented 
and discussed the results in Section IV, while the conclusions are stated in Section V. 




(1) 
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II. The confined Ising ferromagnet with competing fields and the Monte 
Carlo simulation method. 



The Ising model with competing surface fields in a confined geometry of size L x M (L -C M) can be described by 
the following Hamiltonian H: 

M.L M M 

<ij,rnn> i—1 i—1 

where Uy are the Ising spin variables, corresponding to the site of coordinates that may assume two different 

values, namely (Tij = ±1, J > is the coupling constant of the ferromagnet and the first summation of Equation ^ 
runs over all the nearest-neighbor pairs of spins such as 1 < z < Af and 1 < j < L. The second (third) summation 
corresponds to the interaction of the spins placed at the surface layer j = 1 {j = L) of the film where a short range 
surface magnetic field hi {h^) acts. Open boundary conditions are assumed along the M-direction of the film where 
the fields act. 

In this manuscript, only the case of competing surface fields such as hi = —h^, in the absence of any bulk 
magnetic field, is considered. So, hereafter h = \hi \ = will be used generically to specify the short range surface 
magnetic field that is measured in units of the coupling constant J. Under these conditions and for a suitable 
range of fields and temperatures, one may observe the onset of an interface between magnetic domains of opposite 
direction running along the film, as will be discussed in detail below. 

The evolution of the Ising film is simulated using the Glauber dynamics, so a randomly selected spin is fiipped 
with probability p(flip) given by: 

exp(-/3 ■ AH) 
P^^''^^ ^ l + exp(-/3.Ag) 

where AH is the difference between the energy of the would-be new configuration and the old configuration, and 
P = l/ksT is the usual Boltzmann factor. 

The time is measured in Monte Carlo time steps (mcs) , such as during one mcs all L x M spins of the sample 
are fiipped once, in the average. 

The Ising magnet in two dimensions and in absence of any external magnetic field, undergoes a second-order 
order-disorder transition when the temperature is raised from a relatively low initial value. The critical point is the 
so called Onsager critical temperature ksTc/J — 2.269... [36] . 

In order to evaluate the damage according to equation ([l]), first the reference configuration S"^ is generated, 
using the simulation method described above. This can be achieved starting from a random configuration and 
applying the Glauber dynamics during 10'* mcs. Subsequently, a replica of such configuration is created and the 
spins of the central column (i = ^) are fiipped as follows: if the magnetization of the whole sample is greater than 
zero the up spins are fiipped down, otherwise the down spins of the column are fiipped up. Using this procedure the 
perturbed configuration of equation ^ is generated and it is assured that the initial damage (at t = 0) is always 
D{t = 0) < ^ . This kind of perturbation reproduces the effect of a large magnetic field applied at the middle of 
the sample and pointing away to the opposite direction than the actual magnetization of the whole sample. 

The time evolution of the damage D(t), evaluated using equation ([T]), is then followed applying the Glauber 
dynamics to both configurations. At this stage it is essential to use the same random numbers in both systems in 
order to perform the updates [50] . 

Applying the Glauber dynamics to the Ising ferromagnet, it is known that exists a critical temperature for 
damage spreading (To), such as for T > Tjj the damage spreads out over the whole sample while for T < Tjj the 
damage becomes healed after some finite time [50]. Extensive numerical simulations due to Grassberger |5H [52] 
have shown that, in two dimensions and in absence on external magnetic fields, Tn = 0.992Tc. 



III. Discussion of Equilibrium Configurations in Confined Geometries 
with Competing Fields. 

Since the aim of this manuscript is to investigate the role played by the interface between magnetic domains in 
the propagation of damage, it is worth discussing the different equilibrium configurations characteristic of confined 
geometries under different surface field conditions, close to the wetting phase transition. 



3 



The situation originated by the presence of competing surface magnetic fields hi = —h^ (see equation ^) 
can be described in terms of a wetting transition that takes place at a certain -field dependent- critical wetting 
temperature Tyj{h). In fact, for T < T^, a small number of rows parallel to one of the surface of the film have an 
overall magnetization pointing to the same direction than the adjacent surface field. However, the bulk of the film 
has the opposite magnetization (i.e. pointing in the direction of the another competing field). Also the symmetric 
situation is equivalent to the previous one due to the spin-reversal field-reversal symmetry. These configurations, 
where the interface between domains is tightly bounded on the surface, are characteristic of the non-wet phase 
of the system that occurs at low enough temperatures and fields (see figure 1). Increasing the temperature, such 
interface moves farther and farther away from the surface towards the bulk of the film. Just at the interface is 
located in the middle of the film and the system reaches the wet phase for the first time. For T > the interface 
is not longer locaHzed and the system is within the wet phase (see figure 1). 

It should be noticed that the above discussed wetting transition takes place in the thermodynamic limit only 
and the phase diagram has been evaluated exactly by Abraham [37j, yielding 

exp(2J/?) • [cosh(2J/3) - cosh(2/ic/3)] sinh(2J/3), (4) 

where hciT) is the critical surface field (the inverse function of the wetting temperature Tw{h)). 

As discussed above, a well defined wetting transition takes place in the thermodynamic Hmit only. Nevertheless, 
as shown in figure(2) a precursor of this wetting transition can also be observed in confined geometries for finite 
values of L. This precursor is most correctly described in terms of a localization-delocalization transition of the 
interface, which takes place at L— dependent critical temperatures (T„(L, h)). In fact, for T < Tw{L, h) (figure 2(a)) 
one observes coexistence of two phases, each of them having opposite magnetization. Furthermore, the interface 
between domains is located close to one of the surfaces of the film (locaHzed interface; non-wet regime). For 
T > T^{L,h) (figure 2(c)), the wall between domains moves along the i-direction and the system enters to the wet 
regime (delocalized interface). The mean position of the interface remains close to the center of the film just at 
T^{L,h) (figure 2(b)). 

According to the finite-size scaling theory [3l[6] Tyj{L, h) shifts towards Tw{h) according to: 

Tw{oo, h) — Tw{L, h) ^ constant x L~^, (5) 

when the system size L tends to the thermodynamic limit. For further discussions on the locaHzation-delocalization 
transition of the Ising system in confined geometries see e.g. [3, lOj. 

Finally, it is worth discussing the equilibrium configurations obtained in absence of magnetic fields and close 
below to T(7, as shown in figure 2(d). In this case the confined system forms a nonuniform state given by a succession 
of up-spins and down-spins domains, which are essentially ordered at length scales such as the standard correlation 
length {£) remains smaller than the characteristic size of the system (^ < L), but the domain walls are randomly 
placed. The operation of surface magnetic fields favors the parallel spins, and thus suppress this domain formation. 
These facts can be recognized from a visual inspection of the domain patterns shown in figure 2. 

IV. Results and discussion. 

Monte Carlo simulations have been performed using L x M-lattices, for the choices L — 12, 24, 48 and M = 50 x L. 

Preliminary runs indicate that the damage spreading transition can be observed, in principle, for any range of 
temperatures and fields close to the critical wetting line as defined by equation ([4]). However, when the temperature 
is very low the system becomes almost frozen making difficult the acquisition of data with reliable statistics. Also, 
for T > Tc the system is disordered preventing the occurrence of the wetting transition. Taking into account these 
constraints, an extensive study of damage spreading has been performed only close to two particular points of the 
whole wet non-wet phase diagram, namely for \hi\ = = 0.5;T^(ft. — 0.5) = O.SGSTc, and for — 0.90Tc; \hi\ = 
= 0.4291. Additional studies have also been carried out close to = G.75Tc; | = {h^l = 0.6652 and 

= 0.95Tc; \hi\ = ^ 0.3053. 

On view of the obtained results, it is expected that the main conclusions obtained from the study will hold close 
to the whole critical wetting curve. 

In order to characterize the global dynamics of the propagation of the perturbation, the whole damage given 
by the Hamming distance D(t) (see equation ((T|)) and the survival probability P{t), that is the probability that at 
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time t the damage is still propagating, have been evaluated. Also, to study the spatio-temporal evolution of the 
perturbation, the damage profile in the perpendicular (parallel) Px{i,t) {Py{i,t)) direction is defined as follows: 

= ^ E I ^6 ^) - n I , (6) 

and 

M 

where Sf,j[t,T) and S^,j{t,T) are the unperturbed and perturbed equilibrium configurations, respectively. This 
definition of the perpendicular (parallel) damage profile represents the average damaged sites of the i—th -i — 1, .., M 
column { j ~ th - j — 1, .., L row) of the lattice. In order to characterize the profiles it is also useful to determine 
their width {{LUx{t)) and {ujy{t))), as well as their amplitude {{hx{t)) and {hy(t))), respectively. 

The propagation of the damage is studied using two approaches: i) keeping h = 0.5 constant and varying T, 
and ii) keeping T = 0.90 constant and varying h. 



IV. 1 The case h = 0.5, T variable. 

Figure 3 shows log-log plots of the damage D{t) as function of time. It is interesting to note that the locaHzation- 
delocaHzation transition is roughly revealed by the behavior of the damage. In fact, for T < Tw{L), (e.g. for 
T = O.lOTc, O.GbTc and O.GOTc in figure 3(a) ) it is observed that the damage is healed after a relatively short 
time. However for T > Tw{L) (e.g. for T ~ O.SOTc and O.SSTc in figure 3(a)) the spreading of damage is actually 
observed. 

It should be noticed that using finite lattices, the onset of damage spreading is observed at i-dependent "critical" 
temperatures {T£,{L)). As expected, finite-size effects cause the actual transition to become rounded and shifted, 
as e.g. can be observed comparing figures 3(a) and 3(b), which have been obtained using lattices of different size. 
The straight lines observed in figures 3(a) and 3(b) for T > Tu{L) suggest a power-law behavior such as 

D{t) cx (8) 

where -q is an exponent. Disregarding the early time behavior of D{t) (say up to i « 30 mcs), the best fit of the 
data shown in figure 3(a), using equation gives 77 = 0.90 ± 0.02. Simulations performed using bigger lattices 
{L = 24,48, see e.g. figure 3(b)) also give 77 = 0.91 ± 0.01, where in all cases the error bars merely refiect the 
statistical error. 

Notice that damage spreading is also observed for a temperature slightly below T^, (e.g. at T = O.SOTc in figure 
3(a) and T = 0.82Tc in figure 3(b)). This observation is due to the operation of finite size effects causing the 
effective location of the prewetting temperature to be slightly below Tw (see eq.|[5])). Eventually, this result may 
also indicate that the critical temperature for damage spreading Hes below T^, as will be discussed below. 

The survival probability of the damage P{t) is shown in figure 4, as a function of time t. It is observed that 
Pit) tends to vanish for T < Td{L), while it reaches a stationary value for T > Td{L). So, it is possible to define 
the stationary survival probability Pstat as the asymptotic limit Pstat = Pit 00), (see figure 5(a)). In this case, 
finite-size effects are also observed. In fact, the step function expected for the PstatiT) becomes rounded for lattices 
of size L — 12 and M — 601. However, the damage temperature 7z)(L) can be determined as the mean value of the 
rounded step. 

Figure 6(a) shows log-log plots of the width of the damage profile {oJxit)) vs t obtained taking h = 0.5 and 
different temperatures. It is observed that for T > Td(L) ( T < TuiL)) the width exhibits an upward (downward) 
curvature. However, a straight line behavior emerges for a certain temperature that is identified as ToiL) w 0.78. 
Based on these observations, the following scaling law is proposed, 

K(t))cxt", (9) 

where a is an exponent. The best fit of the data gives a — 0.80 ± 0.01. 

It is worth mentioning that a different behavior has been observed for zero surface fields (ft, = 0) close to the 
critical temperature 06], as shown in figure 6(b). In this limit the critical temperature for damage spreading is 
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well known, namely Td(/i = 0, L = oo) ~ 0.992Tc; jsTHH]- However, in this case a crossover between a short-time 
(ST) regime and a long-time (LT) regime is found close to t 350 mcs. For t < 350 mcs the data is consistent 
with asT = 0.517 ± 0.005, while for t > 350 mcs one has ult = 0.79 ± 0.01. This results also suggests that a 
diffusive-like behavior, as expected for a = 1/2, dominates the short-time regime, while the propagation of the 
perturbation becomes faster within the long-time regime. 

The observed difference between the cases h > and h ^ can be understood with the aid of the snapshot 
configurations shown in figures 2(a-c) and 2(d), respectively. In fact, when the interface is most likely localized 
at the center of the strip running parallel to the M-direction, where the surface magnetic fields are applied, one 
expects an uniform propagation of the perturbation as observed in figure (6(a)). However, for h = the occurrence 
of magnetic domains of opposite orientation exhibiting interfaces running perpendicular to the strip determines the 
onset of two different time regimes: on the one hand the short-time regime, characterized by the (slow) propagation 
of the damage inside the bulk of magnetic domains and, on the other hand, the long-time regime with a (faster) 
propagation of the perturbation effectively catalyzed by large fiuctuations occurring at the interfaces of the domains. 

The amplitude of the damage profile in the direction parallel to the interface {hy(t)), as shown in figure 7, also 
exhibits a power-law behavior of the type: 

{hy{t))(Xt^, (10) 

where A is an exponent. A least-square fit of the data (see figure 7) gives A = 0.98 ± 0.01, for T > Tu[L). Also, 
the damage tends to become healed for T < Tu{L). Simulations performed using bigger lattices shows that this 
exponent tends to A w 1 for the asymptotic limit, L oo, e.g. A(L = 24) = 1.03±0.02 and \{L = 48) = 1.04±0.05. 

Figure 8 shows damage profiles Py{i,t) obtained at different times, and close to Td{L). For the absence of 
surface fields and using open boundary conditions (figure 8(a)), the homogeneous propagation of the damage is 
observed. However, considering competing surface fields the damage tends to propagate in the region where the 
interface between magnetic domains is located, namely at the central rows of the film (« ~ L/2), as shown in figure 
8(b). This fact is in qualitative agreement with the results discussed for {hy(t)). Moreover, it is observed that the 
effect of the magnetic fields at the walls is to slow down the propagation of the damage. These observations are 
consistent with the fact that large fluctuations capable to catalyze the propagation of the damage are present along 
the interface, while at the walls and due to the appHed magnetic flelds, the orientation of the spins is more stable. 

It should be noted that due to the finite size of the lattice, the width of the damage profile (wy(i)), remains 
constant close to -j, preventing a quantitative analysis. 

Based on the study of D{t), P{t) and the damage profiles, the location of the L— dependent damage spreading 
critical points Td{L), for a fixed field h = 0.5, have been estimated as follows: Td{L ~ 12) — 0.79(l)Tc, Td(I/ = 
24) = 0.81(l)rc and Td{L = 48) = 0.83(l)Tc. Based on these results and assuming that Td{L) obeys a scaHng 
law such as equation |[5]), the extrapolation Td{oo) — 0.84 ± 0.01 has been obtained, as shown in figure 9. This 
value is close (but smaller) to the wetting critical temperature Tw{h — 0.5) — 0.863, as evaluated using eq.|[4]) (see 
also figure 1). It is known that one method to locate the L— dependent critical temperatures for the localization- 
delocaHzation transition (Tu,(L)), which is the precursor of the true wetting temperature in the thermodynamic 
limit, is the evaluation of the maximum of the susceptibility, i.e. the fluctuations of the order parameter [3]. Figure 
9 also shows that a plot of Tw{L) vs (see equation ([5])) extrapolates to Tw{L oo) ~ 0.866, in good agreement 
with the exact value. Furthermore, one has that T]j{L) < Tw{L). 

At this stage it is concluded that, on the one hand, it is hard to establish an unambiguous difference between 
Td and T^, while on the other hand, it is out of any doubt that the onset of damage propagation is related to 
the locaHzation-delocalization transition that undergoes the interface between competing domains, so that Td is 
effectively shifted to temperatures far below from Tc- 

IV.2 The case T = 0.9Tc, and h variable. 

In this case, the obtained results are fully consistent with those obtained keeping the field constant and varying 
the temperature. In fact, the total damage D{t) grows up for h > hr){L) with exponent rj = 0.91 ± 0.02 (see figure 
10), where hu{L) is L— dependent critical field for damage spreading. Also, it is observed that for lattices of size 
L = 12, the damage spreading transition is markedly rounded by finite-size effects (figure 10(a)), while for L = 24 
that transition is better defined (figure 10(b)). 

The stationary survival probability of the damage Pstat can also be determined, as shown in figure 5(b) as a 
function of the surface field h. The fact that the stepped shape of Pstat is better defined in figure 5(b) than in figure 
5(a) is an intriguing feature that remains to be clarified. 
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The width of the damage profile also behaves according to a power law (see equation ([8]) with an exponent 
a = 0.81 ± 0.01 (see figure 11), in full agreement with our previous estimation given by a = 0.80 ± 0.01. 

It is found that the amplitude of the profile, {hy{t)), also exhibits a power law close to (figure 12). In fact, the 
exponent A (see equation ifTOj) . has been determined given: A = 1.02 ± 0.01, in full agreement with our estimation 
obtained keeping constant the surface fields, namely A = 0.98 ± 0.01. 

It has been observed that the amplitude along the parallel direction and the width of the profiles across the 
perpendicular direction {{hx(t)) and {ujy{t)), respectively), exhibit the same behavior than in the previous case with 
constant surface fields, so the corresponding figures are omitted for the sake of space. 

Finally, our estimations of the size dependent damage spreading critical fields hD(L) are given by: /i_d(L = 
12) = 0.31(2), hD[L = 24) = 0.34(1) and hoiL = 48) = 0.35(1). 

By comparison to equation ^ the following finite size scaHng anzats can be expected to hold for the wetting 
critical field 

hc{L) — hc{oo) ^ constant x L~^, (11) 

where ftc(oo) is the position of the critical point in the thermodynamic limit. Plotting the obtained data for 
the damage spreading critical fields according to equation (flTI) (figure 13), the extrapolation to L ^ oo gives 
hoiL — oo) — 0.37 ± 0.02, while the exact value for the wetting critical point at T = 0.9Tc, as obtained using 
equation |(4]), is ft-c = 0.42906. As in the previous case, the location of the damage spreading critical point lies 
slightly inside the non-wet region of the phase diagram (see figure 1). 

IV. 2 Complementary studies. 

In order to confirm the findings of the previous sections, additional simulations were performed at T = 0.75Tc 
and T = 0.95Tc and varying the surface fields. All the obtained results are consistent with the previous findings. 
In summary, applying the scaling anzats given by Equation (fTTj) . the critical field for damage spreading /i£)(oo) = 
0.62±0.01 (ft-D(oo) — 0.15±0.01) was obtained, as compared with the exact result for the wetting critical point (e.g. 
equation ([4])) given by hdT = 0.75Tc) = 0.6652 {hc{T = 0.95Tc) = 0.3053). Notice that the former extrapolation 
is shown in figure 13. 

V. Conclusions 

Based on an extensive numerical study of the damage propagation close to the wetting transition of the Ising model 
with competing fields it is possible to draw the following main conclusions: 

i) The exponent governing the spreading of the total damage {D{t) cx f, see Eq.®), in the presence of interfaces 
in the direction parallel to the propagation of the perturbation, namely rj = 0.90 ±0.02, is greater than the obtained 
in absence of interface, namely r] = 0.471 ± 0.005 [HI |46]. So, clearly the interface causes the enhancement of the 
propagation of the perturbation. 

From a more qualitative point of view, it is expected that the propagation of damage may be favored by interfaces 
because precisely around these regions of the sample one has the largest fiuctuations in the orientation of the spins. 
Within this scenario the spreading of the damage is expected to be faster along the direction parallel to the interface 
where the fiuctuations can uniformly be propagated. In contrast, the spreading along the direction perpendicular 
to the interface may be slowed down by the periods where the damage has to cross over the bulk of the domains 

m- 

ii) Clear evidence on the shift of the critical damage temperature is reported. In fact, while it is well known that 
Toih = 0) « 0.992Tc in the absence of surface magnetic fields [5ll[52], our finding show that Toih > Q) < Toih = 0) 
when competing surface fields are appHed. Furthermore, all the damage spreading critical points extrapolated to 
the thermodynamic limit He systematically within the non-wet phase of the wetting phase diagram (figure 1). These 
results suggest that the wetting and the damage spreading critical curves are different and point out that the former 
should be slightly shifted towards the non-wet phase. 

We expect that this numerical work may contributed to the understanding of the role of interfaces in the 
propagation of damage. The overall understanding of these complex physical phenomena addresses a significant 
theoretical and experimental challenge. Such kind of studies are stimulated by many practical and technological 
applications as e.g. the design and construction of high quality magnetic materials. 
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Figure 1: Wetting phase diagram {h vs T/Tc) corresponding to a semi-infinite Ising strip with opposite short range 
surface fields. The solid curve is the exact critical line, as obtained by Abraham [37] . The vertical dashed line is the 
boundary between the disordered phase at right hand side and the ordered one at the left hand side. The results 
obtained for the damage spreading transition are shown as full circles, and the dotted line has been drawn to guide 
the eyes. Crosses indicate the three points, of the phase diagram, where the snapshots pictures shown in figure 2 
were obtained. 
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Figure 2: Snapshots configurations obtained after 10^ mcs in confined geometries of size L = 24 and M = 1200. 
(a)-(c) Typical spin configurations when short-range fields with opposite signs are appHed at bottom and top files 
of the lattice, respectively. In these cases, the snapshots are taken at T = 0.80 Tp and difi"erent surface fields: (a) 
h = 0.4, within the non-wet phase; (b) h — 0.6, near the critical wetting curve, and (c) h — 0.8, within the wet 
phase. The location of these points are shown in the phase diagram of figure 1. (d) Typical equilibirum configuration 
obtained in absence of magnetic fields, applying open boundary conditions, and close to Tc, T = O.OOTc. Note 
that the horizontal coordinate has been reduced by a factor of five in comparison with the vertical one for the sake 
of clarity of the picture. Sites taken by downs spins are shown in black while up spins are left white. 
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Figure 3: Log-log plots of D{t) vs t obtained at different temperatures as indicated in the figure, applying short- 
range fields of magnitude h = 0.5, and using lattices of size: (a) L = 12, M = 601, (b) L = 24,M = 1201. The 
dashed lines have slopes t] = 0.90 and have been drawn for the sake of comparison. 
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Figure 4: Log-log plots of P{t) vs t obtained at different temperatures as indicated in the figure, and using lattices 
of size L = 12 and M = 601. Short-range fields of magnitude h = 0.5 are applied at the walls of the sample. 
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Figure 5: Plots of the stationary survival probability Pgtat = oo) as obtained using lattices of size L = 12, 

Al = 601. (a) PsLai vs T/Tc obtained taking short-range fields of magnitude h = 0.5 (the dotted line has been drawn 
to guide the eyes). The vertical dashed line indicates the exact critical wetting point: Tw{h^ = 0.5) = 0.863Tc. 
(b) Pstat h obtained taking a fixed temperature T = O.QOTcr (the dotted line has been drawn to guide the eyes). 
The vertical dashed line corresponds to the exact wetting^critical point: /i«,(T^ = O.QOTc) = 0.4291. 




Figure 6: Log-log plots of < ujx{t) > vst obtained at different temperatures, as indicated in the figure, using lattices 
of size L — 12 M = 601. (a) Results obtained taking short-range fields of magnitude h = 0.5. The dashed line has 
slope a — 0.80. (b) Results obtained in the absence of magnetic fields [46] and using lattices of different size, with 
M = 50 X L, as indicated in the figure. The dashed line has slope asT = 0.52, for t < 350, and the full line has 
slope olt — 0.78, for t > 350, respectively. More detailSj^i^ the text. 
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Figure 7: Log-log plots of < hy{t) > vst obtained at different temperatures, as indicated in the figure, using lattices 
of size L = 12 M = 601, and applying short-range fields of magnitude h = 0.5. The dashed line has slope A = 0.98. 
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Figure 8: Plots of the damage profile measured along the L— direction and obtained at different times as indicated 
in the figure. Results are averaged over 10^ samples, using lattices of size L = 12, M = 601. (a) T = 0.98Tc and 
open boundary conditions, (b) T = 0.861Tc and short-range fields \hi \ = = 0.5. 
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Figure 9: Damage spreading critical points Td{L) as a function of L^^ , obtained for lattices of different size 
(L — 12, 24, 48, with M = 50 x L), and applying short range surfaces fields of magnitude h = 0.5 (full circles). The 
dashed line extrapolates to the critical temperature T^jTc — 0.84 ± 0.01. Full squares correspond to the interface 
localization-delocalization transitions T^(V) as obtained for h — 0.5 and using the criteria of the maximum of the 
susceptibility [5]. The full line extrapolates to the wetting critical temperature T^jTc = 0.866. 
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Figure 10: Log-log plots of -D(t) vs t obtained using short-range magnetic fields of different magnitude as indicated 
in the figure, at fixed temperature T = 0.90Tc, and using lattices of size: (a) L = 12, M = 601, (b) L = 24, 
M = 1201. The dashed lines have slopes 77 ~ 0.90. 
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Figure 11: Log- log plots of < U!x{t) > vs t obtained for short-range magnetic fields of different magnitude, as 
indicated in the figure, using lattices of size L = 12 M = 601, and at fixed temperature T = O.QOTc. The dashed 
line has slope a — 0.81. 
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Figure 12: Log- log plots of < hyif) > vs t obtained at different short-range magnetic fields, as indicated in the 
figure, using lattices of size L = 12 M = 601, and at fixed temperature T = 0.90Tc. The dashed line has slope 
A = 1.02. 
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Figure 13: Damage spreading critical fields hoiL) as a function of L ^, obtained for lattice of different size 
{L — 12, 24, 48), and at different temperatures, as indicated in the figure. The dashed (full) line extrapolates to the 
critical field ^ oo,T = 0.75Tc) = 0.62 ± 0.01 ^ oo,T = 0.90Tc) = 0.37 ±0.02). Both lines have the 

same slope. 
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